Spin-excitations of the quantum Hall ferromagnet of composite fermions 
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The spin-excitations of a fractional quantum Hall system are evaluated within a bosonization 
approach. In a first step, we generalize Murthy and Shankar's Hamiltonian theory of the fractional 
quantum Hall effect to the case of composite fermions with an extra discrete degree of freedom. 
Here, we mainly investigate the spin degrees of freedom, but the proposed formalism may be useful 
also in the study of bilayer quantum-Hall systems, where the layer index may formally be treated 
as an isospin. In a second step, we apply a bosonization scheme, recently developed for the study 
of the two-dimensional electron gas, to the interacting composite-fermion Hamiltonian. The disper- 
sion of the bosons, which represent quasiparticle-quasihole excitations, is analytically evaluated for 
fractional quantum Hall systems at i/ = 1/3 and u = 1/5. The finite width of the two-dimensional 
electron gas is also taken into account explicitly. Furthermore, we consider the interacting bosonic 
model and calculate the lowest-energy state for two bosons. In addition to a continuum describing 
scattering states, we find a bound-state of two bosons. This state is interpreted as a pair excitation, 
which consists of a skyrmion of composite fermions and an antiskyrmion of composite fermions. The 
dispersion relation of the two-boson state is evaluated for u = 1/3 and v — 1/5. Finally, we show 
that our theory provides the microscopic basis for a phenomenological non-linear sigma-model for 
studying the skyrmion of composite fermions. 

PACS numbers; 71.10.Pm, 71.70.Di, 73.43.Cd, 73.43.Lp 
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I. INTRODUCTION 

Quantum Hall physics - the study of two-dimensional 
(2D) electrons in a strong magnetic field - has revealed a 
lot of unexpected phenomena during the last 25 years.^ 
Apart from the integer and fractional quantum Hall ef- 
fects (IQHE and FQHE, respectively), which have been 
identified as macroscopic quantum phenomena, exotic 
topological spin excitations are displayed by these sys- 
tems. The physical properties of quantum Hall systems 
are governed by the Landau quantization; the energy of 
2D electrons with a band mass and charge — e in a 
perpendicular magnetic field B is quantized into equidis- 
tant energy levels, the so-called Landau levels (LLs), with 
a level separation ?iz/;c = TieB/mBC. The large LL degen- 
eracy is characterized by the fiux density = l/2iTlg, 
given in terms of the magnetic length I b — ^/fic/eB, and 
the LL filling is thus defined as the ratio v — riei/n^ of 
the electronic density rie; and n^. Due to the Zeeman 
effect, each LL n is split into two spin-branches with 
an energy separation g* ^bB, where g* is the effective 
Lande factor (g* — —0.44 for GaAs), and hb = 2he/m is 
the Bohr magneton. Because the band mass is reduced 
with respect to the bare electron mass m (ms = 0.068m 
for GaAs), the spin-branch separation is about 70 times 
smaller than the LL separation. 

The IQHE may be understood in a one-particle pic- 
ture, in which the Coulomb interaction is only a small 



perturbation. If = A^, with integral A^, the ground 
state is non-degenerate and one has to provide a finite 
energy to promote an electron to the next higher level. 
Due to the localization of additional electrons by residual 
impurities in the sample, the Hall resistance remains at 
its quantized value Rh = hje^N over a certain range of 
the magnetic field. This plateau in the Hall resistance is 
accompanied by a vanishing longitudinal resistance, and 
both are the signature of the quantum Hall effect. If the 
lower spin branch of the n-th LL is completely filled, one 
has A^ = 2n -f 1, and A^ = 2n -I- 2 in the case of complete 
filling of both spin branches. 

The one-particle picture ceases to be valid at par- 
tial filling V ^ N , where one is confronted with the 
LL degeneracy. In this limit, electrons in a partially 
filled level are strongly correlated due to their mutual 
Coulomb repulsion, which constitutes the relevant en- 
ergy scale. In the lowest LL, the formation of composite 
fermions (CPs) leads to incompressible quantum liquids 
aiv = p/(2sp-|- 1).2>3 The PQHE, displayed by these liq- 
uids, may be interpreted as an IQHE of CPs. The physi- 
cal properties of CPs have been extensively studied in the 
framework of Jain's wave- function approach, ^■'^ which is a 
generalization of Laughlin's trial wave functions.^ A com- 
plementary field-theoretical formalism, based on Chern- 
Simons transformations, has been proposed by Lopez and 
Pradkin.^ The latter has been particularly successful in 
the description of the metallic state at v = 1/2, which 
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may be interpreted as a Fermi sea of CFs.^ More recently, 
Murthy and Shankar have developed a Hamiltonian the- 
ory of the FQHE, a second-quantized approachii It sys- 
tematically accounts for the mechanism of how these 
excitations are formed and what are their constituents, 
namely, electrons or holes bound to an even number of 
vortex-like objects. 

Although the main features of the IQHE may be un- 
derstood in the framework of a noninteracting model of 
spin-polarized electrons, the study of spin-excitations at 
1/ = 1, when only one spin-branch of the lowest LL 
(n — 0) is completely filled, requires a more complete 
treatment of the problem. Indeed, the lowest-energy ex- 
citations involve a spin reversal, and the spin degree of 
freedom must be considered explicitly. In addition to 
spin-wave excitations, which are gapped due to the Zee- 
man effect, one finds topological skyrmion excitations)^ 
if the Coulomb interaction is large with respect to the 
Zeeman gap, it is energetically favorable to distribute 
the spin reversal over a group of neighboring spinsAiSiii 
These skyrmion excitations give rise to unusual spin 
polarizations.^^ More recently, the different spin exci- 
tations at v — 1 have been investigated in a bosoniza- 
tion approach. Although the bosonization method of 
fermionic systems is well established for one-dimensional 
electron systemsji^^ only few generalizations to higher 
dimensions have been proposed j^^i^^ One particular ex- 
ample is a recently developed bosonization scheme for 
the 2D electron gas (2DEG) in a magnetic fieldiiii£ The 
extension presented in Ref. Il^i in particular, allows for 
the treatment of complex many-body structures, such 
as small skyrmion-antiskyrmion pairs in terms of bound 
states of the bosonic excitations. 

Here, we generalize the Hamiltonian theory of the 
FQHE of Ref. to incorporate spin or any other discrete 
degree of freedom. Although in this paper we apply our 
theory to the investigation of spin excitations, the devel- 
oped formalism has a wider range of application and may 
also be useful for further studies on bilayer quantum Hall 
systems, for which the layer index is formally treated as 
an isospinpiiii For certain values of the magnetic field, 
the ground state may be viewed as a = 1 state of 
CFs, and it is therefore natural to test the bosonization 
scheme for these cases. The method allows us to inves- 
tigate the CF excitation spectrum and to describe, for 
example, skyrmion-like excitations at filling factors such 
as v = 1/3. Indeed, optically pumped nuclear magnetic 
resonance measurements, which give information about 
the spin polarization of the 2DEG, indicate the existence 
of such excitation around ly = 1/3."'^^ 

The outline of the paper is the following. The general- 
ization of the Hamiltonian theory to the case of fermions 
with an extra discrete degree of freedom is presented 
in Sec. II. In Sec. Ill, we apply the previously devel- 
oped bosonization method^^ to the obtained CF Hamil- 
tonian and study the dispersion relation of the resulting 
bosons. Bound states of two of these elementary excita- 
tions, which we interpret as small skyrmion-antiskyrmion 



pairs, are also considered. Sec. IV is devoted to the semi- 
classical limit of our approach, from which we show that 
an effective Lagrangian, similar to the one suggested by 
Sondhi et al.^ accounts for the long- wavelength prop- 
erties of the system. A summary of the main results is 
finally presented in Sec. V. 



II. HAMILTONIAN THEORY OF CFS WITH 
SPIN 

Our first aim is to derive the Hamiltonian theory for 
electrons with spin in the lowest LL and construct the 
CF basis. In principle, the generahzation of the alge- 
braic properties of the Murthy and Shankar model [Ref. 
01 to the case of discrete degrees of freedom, such as the 
spin or the isospin in the context of bilayer quantum Hall 
systems, may be written down intuitively. In this section, 
however, they are derived, with the help of a set of trans- 
formations, from the microscopic model of an interacting 
2DEG, with spin a =t (+1) or | (-1) in a perpendicular 
magnetic field Bgz = V x A. The Hamiltonian of the 
system reads, in terms of the electron fields V'a-Cr) and 

V4(r) {n = c = i) 

H = Ho + Hz + Hi, (1) 
where the orbital part is 

^^ = ^11 j d'r^Ur) + eA{r)f V^,(r), (2) 

the Zeeman term is 

Hz = \9*liBBY,cj I d2rV4(r)^.(r), (3) 

and the interaction Hamiltonian is 

<T.(T' 

(4) 

which couples electrons with different spin orientation. 
Here, V{y) = e^/er is the isotropic Coulomb interaction 
potential, with r — |r|, and e is the dielectric constant of 
the host semiconductor. 



A. Microscopic Theory 

Following closely the steps of Ref. 0, we perform a 
Chern-Simons transformation on the electron fields, de- 
fined by 

; / \ -2js„ r <fr' e(T-r')p„(r') iCSt \ /rr\ 

ip„{v)^e 'J ^ 'ip^^{r), (5) 

where 6{r) is the angle between the vector r and the x- 
axis. The density operator reads /Oo-(r) = ■(/'i(r)'0cr(r) = 
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tp^^^{r)ip^^{r). Thus, the Zeeman Hz and the inter- 
action Hj terms in Eq. Q remain invariant under the 
above transformation. In order to describe fermionic 
fields, which is the case of interest here, must be in- 
teger. Note that this transformation breaks the SU{2) 
symmetry exphcitly unless S| = Sj. For the symmet- 
ric case, a Hamiltonian theory for bilayer systems has 
very recently been established by Stanic and Milovanovic 
who investigated composite bosons at a total filling fac- 
tor 1/ — 1.^" Here, we discuss the more general case in 
which the electronic density of "f-electrons is not neces- 
sarily the same as that of |-electrons. The transforma- 
tion generates the two-component spinor Chern-Simons 
vector potential a^'^(r) with 



X a^^» 



-2so-(/>o/Ocr(i")e2 



(6) 



in terms of the flux quantum (j>o = 27r/e. 

In order to decouple the electronic degrees of freedom 
and the complicated fluctuations of the Chern-Simons 
field, Murthy and Shankar introduced an auxiliary vec- 
tor potential^ a(j(q), which was inspired by the plasmon 
description proposed by Bohm and Pines for the elec- 
tron gas.^^ If the spin degree of freedom is taken into 
account, this vector potential consists of two spinor com- 
ponents, both of which are chosen transverse to satisfy 
the Coulomb gauge, V • aa.(r) = 0. For its Fourier com- 
ponents one thus finds ao-(q) = —ie^ x eqa^{q), where 
Bq = q/|q| is the unit vector in the direction of propaga- 
tion of the field. Its conjugate field is chosen longitudinal, 
Po-(q) = i^qPail), with the commutation relations 



r(q),P<T'(-q')] = iSa,a'Sq,q'. 



(7) 



The Hilbert space is enlarged by these non-physical de- 
grees of freedom, and one therefore has to impose the 
constraint ao-(q)|0) = on the physical states The 
orbital part Q of the Hamiltonian thus becomes 



X [n^ + e:a^^(r) :+ea^(r)]'^^^(r), (8) 

where the mean-field part of the Chern-Simons vector po- 
tential (a^'^(r)) has been absorbed into a generalized mo- 
mentum Tier = -iV+e [A(r) + (a^'^(r))] , and : a^'^(r) : 
denotes its normal-ordered fluctuations. The conjugate 
field Pa (q) generates the translations of the vector poten- 
tial, and one may therefore choose the unitary transfor- 
mation 'ip^^{r) = Uatpa^ {ic), in terms of the composite- 
particle field ipa^{r), with 



Hamiltonian (jSJ without the short-range terms related to 
: a^^{q) : thus reads 



(10) 



with the Hamiltonian density 

n! 



^ ^a„(r)+4^s„P<,(r)]' 



2m B 2m b 
+— n,, • [a^(r) + 47rs^P^(r) 

TUB 



(11) 



The second term is the Hamiltonian density of a quantum 
mechanical harmonic oscillator, as may be seen from the 
commutation relations (TJ. In the absence of the third 
term of Eq. , which couples the electric current den- 
sity and the harmonic oscillators, one may write down 
the field as a product of the ground-state wavefunction 
of the oscillators x(r) and 4>*{r), where (f>*{r) is the N- 
particle wavefunction of charged particles in an effective 
magnetic field 



Vx [A(r) + (a^^(r))] =SX, 



(12) 



with B* = B - 2sa4>on'^i = Ba/{2saPa + 1), in terms of 
the average density n^j of electrons with spin a. In the 
absence of correlations between electrons with different 
spin, the filling of each spin branch may thus be char- 
acterized by a renormalized filling factor, z/* = ng;27rZ*^, 
with the new magnetic length /* = 1/ ^JeB*. The oscil- 
lator wavefunction becomes, with the help of the trans- 
formed constraint, [qa„{q) — ^■KS„pc,{q)] \(f)) = 



X(r) = exp 



V 27rs<,^/9<^(-q)— (5p<,(q) 



The exponent may be rewritten in terms of the complex 



coordinates Zj 



■ iyj for spin-up and = -t- iy^ 



for spin-down particles^^ 

2nSa5pa{~q)^Spa{q} 



CT,q 



2s| In \zi 



where — 2sa-Pa / {2saPa + 1), at integer values of v* = 
Pa- The complete wavefunction therefore becomes 



Ua =exp i^-PCT(-q)^^(5/3(q) , (9) '^Piz^,Wf,) = Y[{zi - Zjf"^ Y[ {w,, - w^f'^ (t)*{z,,Wf,), 



to make the long-range {q < Q) fluctuations of the 
Chern-Simons vector potential vanish. The transformed 



(13) 

where the Gaussian factors have been absorbed into the 
electronic part of the wavefunction (/)*(zi,w^). If there 
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are no correlations between particles with different spin 
orientation, i.e. in the absence of interactions, the elec- 
tronic part of the wavefunction may be written as a prod- 
uct, (j)*izi,Wfj_) = 0l,^(zi)0^^(w^), and Eq. (H^Jl thus rep- 
resents a product of two Laughlin wavefunctions^ (for 
Per — 1) or unprojected Jain wavefunctions^ (for p^^ > 1) 
for the two spin orientations. If one includes correlations 
between particles with different spin orientation, one may 
use the ansatz 

with integral m. For = 1, one recovers Halperin's 
wavefunctions)^ if 



and similarly 



w^) e 



In this case the filling factors of the two spin branches he- 
come i^cr — (2so- — m)/(s|S|— m^). For larger values of , 
the wavefunction contains components in higher LLs 
and is thus no longer analytic, as required by the lowest 
LL condition.^'* Eq. must therefore be projected into 
the lowest LL, and the final trial wavefunctions become 



i<j nKf 

xYl{zi-w^)"'(l)l^{z,)(j)l^{Wf,), (14) 

where V denotes the projection to the lowest LL. They 
may be interpreted as a CF generalization of Halperin's 
wavefunctions.^^ 

Up to this point, the coupling term between the elec- 
tronic and oscillator degrees of freedom in the Hamil- 
tonian density pil) has been neglected. Murthy and 
Shankar have proposed a decoupling procedure, which 
also applies to the case of electrons with spin, with the 
only difference that two unitary transformations Ua- are 
required for the two spin orientations. The reader is re- 
ferred to their review (Ref. 0) for the details of the de- 
coupling transformations, and we will limit ourselves to 
the presentation of the results. With the choice of equal 
number of oscillators and particles for each spin orienta- 
tion, the kinetic term in the Hamiltonian vanishes, and 
the model in the limit of small wave vectors becomes 

CT.q a 

-^Xm^o'^'^)^'^^^^)^'^'^^) (15) 

q a. a' 



where A^{q) = [ao-(q) -I- AmSaPaifl)] /y/SnSa, 7 
g* ^bB is the Zeeman energy, and 



"0 



(16) 



is the interaction potential in the lowest LL. The trans- 
formed density operators may be expressed in first quan- 
tization with the help of the position of the j-th particle 
Tj and its generalized momentum 11^, 



p.(q) = ^e-''i"-^ 



1 



(17) 

where terms of order 0{q^) have been omitted, and the 
oscillators are chosen to remain in their ground state, 
(Acr(q)) — (Aj.(q)) = 0. The transformed constraint is 
given by 



with 
X.(q)=^e-^'i- 



1 + 



""0,(1 + 0,) 



(18) 



(19) 



To construct the theory at larger wave vectors, Shankar 
had the ingenious insight that if one considers the terms 
in the expressions H17|) and H19|) as the first terms of a 
series expansion of an exponential, 



5,(q)=5:e--<' 



and 



X 



.(q)=5:e--<^. 



both operators find a compelling physical 
interpretation.^'"^ The operator 



Rl^) = r - 



{l + Ca) 



may be interpreted as the electronic guiding center oper- 
ator, and its components satisfy the correct commutation 



relations. 



j^^e) y{e) 



-il%Scr ct', where the number in- 



dex has been omitted. This leads to the algebra for the 
Fourier components of the density in the lowest LIjS& 



[/9,(q),p,'(q')] = Sa,a'2isin 



q A q'/^ 



p<T(q + q')> 



(20) 

where we have defined q A q' = (q x q')z- In the same 
manner, the components of the operator 



Til^^ = r - 



/2 



-e^ X n„ 



satisfy 
[x<T(q),X^'(q')] = 



c,(l- 

i^s'^cr,,'/'^^' which leads to 



qAq'; 



2c? 



^x.(q + q') 



(21) 
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The operator Xcri^): which characterizes the constraint, 
may thus be interpreted as the density operator of a sec- 
ond type of particle with charge — c^, in units of the 
electronic charge. However, this particle, the so-called 
pseudo-vortex, only exists in the enlarged Hilbert space 
because it is an excitation of the correlated electron liq- 
uid and must be described in terms of the electronic de- 
grees of freedom. A similar algebra has been investigated 
by Pasquier and Haldane in the description of bosons at 
1/ = Im^ The final model Hamiltonian thus reads 

q a, a' 

with the constraint H18|l and the algebras H2U|) and H21|l 
for the electronic and the pseudo-vortex density, respec- 
tively. In the case of bilayer systems, the interaction 
between particles in the same layer is different (stronger) 
than that between particles in the different layers. This 
breaks the SU (2) symmetry, and one thus has to replace 



vo{q) 



(<?)- 



B. Model in the CF-Basis 



In order to construct the CF representation of the 
model, Murthy and Shankar proposed the "preferred" 
combination of electronic and pseudo-vortex density^*2^ 



(q), 



(23) 



which plays the role of the CF density. This approx- 
imation allows one to omit the constraint in the case 
of a gapped ground statCj^ whereas it fails in the limit 
v = 1/2, where the system becomes compressible® and 
where the constraint has to be taken into account ex- 
plicitly in a conserving approximationiSS The CF basis is 
introduced after a variable transformation of the guiding- 
center coordinates R^^^ and Rct^\ 

R(.^) - R, + rj^c, R^") = R. + Va/ca. 

The new variables play the role of the CF guiding center 
(Rcr) and the CF cyclotron variable (ry^.) and commute 
with each other. Their components, however, satisfy the 
commutation relations 



in terms of the CF magnetic length /* = Ib/ y/T~c^. 
The CF density operator therefore becomes 



n,n' ;m,m' 



where cj^ ^.^ creates a CF with spin a in the state |n, m). 
The quantum number n denotes the CF-LL, and m indi- 
cates the CF guiding-center state. The matrix elements 
are given by^^ 



G„ 



( ]*\ - ( -i{qx + iqy)l*a ^ 

'^'^ "^ ~ \ ml[ V2 



(25) 



for m > m' , and 



(26) 



2(l-n+n') -{qlB?/2clTn-n' ( <f}^\ 



for n > n' . 

We now restrict the Hilbert space to the lowest CF- 
LL. In this case, the CF density operator p^{r) reads 
(pf (q))„=o = F.^IniQ)pA<l), where 

(27) 

is the CF form factor, p| = p| = 1, S| = S| = s, 
and therefore C| = = c and I* = I* = I* (symmet- 
ric model). The electronic filling factor is thus given by 
h' ~ l/(2s + 1), and the projected CF density operator is 



Pain) = X! ^rn,ni'{ql*)cln^aCm',a, 



(28) 



where cj„ ^ creates a spin a composite fermion in the 
lowest CF-LL and where we have omitted the level in- 
dex. The projected CF density operators satisfy a similar 
algebra as the projected electron density operators^SiS 



[p<T(q):P^'(q') 



<5,,,,2zsinf3^^^)^,(q + q') 



(29) 

just with the magnetic length Is replaced by the CF mag- 
netic length I*. The Hamiltonian, which describes the 
low-energy excitations of the FQHE system at i/* = 1 is 
thus given by 



H 



a, a' q 

Hz + Hi 



(30) 

Note the similarity with the original Hamiltonian H22|l : 
(24) the product VQ{q) \F^£Q{q)\ plays the role of an effective 



6 




n = 



FIG. 1: Schematic representation of the ground state of the 
2DEG at = 1 (the quantum Hall ferromagnet of CFs). 



CF interaction potential. In this sense, the FQHE system 
at V = l/(2s + 1) may be treated exactly in the same 
manner as the 2DEG at v = 1. 

In principle, one may also describe higher CF-LLs 
with only one filled spin branch. If the completely filled 
lower levels are assumed to form a homogeneous inert 
background, the low-energy excitations of the system at 
1/* = 2n+ 1 are described by the same Hamiltonian H3U|) 
now in terms of the CF form factor of the n-th level, 
Pn^il) = Pn,n{q)- Here, however, we concentrate only 
on the case n = 0, and we neglect the presence of higher 
CF-LLs. This assumption becomes valid if the spin split- 
ting is smaller than the typical CF-LL splitting (formally, 
we consider the limit 7^0). 



III. BOSONIZATION OF THE MODEL 

As the model Hamiltonian (|30|l for the FQHE system 
at = l/{2s + 1) resembles the one of the 2DEG at 
v ^ 1, it will be studied with the aid of the bosonization 
method for the 2DEG at v = 1 recently developed by 
two of us^ More precisely, within this bosonization ap- 
proach, we can study the elementary neutral excitations 
(spin-waves, quasiparticle-quasihole pair) of the system. 
Here, we will quote the main results in order to apply the 
bosonization method to the case of CFs at i^* = 1. For 
more details, see Ref. 

In the same way as it was done for the 2DEG at v = 1, 
we first consider the noninteracting term of the Hamil- 
tonian (|30|l . namely the Zeeman term restricted to the 
lowest CF-LL, in order to define the bosonic operators. 
Kt V* = 1, the ground state of the noninteracting system 
is the quantum Hall ferromagnet of CFs, jFAf*^^), which 
is illustrated in Fig. ^ and which can be written as 

n;-i 

\FM^^)= n 4.t|0)' (31) 

m=0 

where |0) is the vacuum state and is the degeneracy 
of each CF-LL. The neutral elementary excitation of the 
system corresponds to a spin reversal within the lowest 
CF-LL as one spin up CF is destroyed and one spin down 
CF is created in the system. This excited state may be 
built by applying the spin density operator S~ = Sx—iSy 
to the ground state 



We define the projected CF spin density operators in 
the same manner as for the projected CF density opera- 
tors [Eq. PI)]. 

-^q = Gm,m'(9^*)Cm |C„i/ I (32) 

5q = ^ G™,„'(gr)c|„^c™.| (33) 

'm,m' 

The commutation relation between and may thus 
be written in terms of CF density operators, 

[§q , S^.] = e-^^'^^''"'%{ci + q')-e-'''^'»''"/'px(q + ^O- 

Although the above commutation relation is quite dif- 
ferent from the canonical commutation relation between 
bosonic operators, its average value in the ground state of 
the noninteracting system [Eq. (|31|l ] is not. Notice that, 
the average values of (q) and (q) in the uniform state 
ED are (Pt(q + q')) = ^^W.o and (^j,(q + q')) = 0, 
respectively, where — A/2ttI*'^ is the number of states 
in one CF-LL, in terms of the total surface A. Therefore, 
the average value of the above commutator is 

([5q,5;,])=iV;<5q+q,^0. (34) 

With the help of the expression H34|) , we can define the 
following bosonic operators 

6q = {n;)-''^~s\, (35) 

h\ - {N;)-"^SI. (36) 

Here, the commutation relation between the operators 5q 
and is approximated by its average value in the ground 
state (|31|) . i.e. 

[V^^']«([^,^^'])='5q.q', (37) 

in analogy with Tomonaga's approach to the one- 
dimensional electron gas.^° This is the main approxima- 
tion of the bosonization method. 

The quantum Hall ferromagnet of CFs \FM'-^^) is the 
boson vacuum as the action of Cm' x on this state is 
equal to zero. Therefore, the eigenvectors of the bosonic 
Hilbert space are 

\{n^})^Y['-^\FM^^), (38) 
q V q 

with riq > and '^q < ■ addition, note that the 
state b'l^FM'-^^) is a linear combination of particle- hole 
excitations. 

One also finds the bosonic representation of the pro- 
jected CF density operators p^[q)^ 

pM) = ^..T^q.oA^; - a^e— '^^'^''^/^fo^^^fok. (39) 

k 
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One can see that the CF density operator is quadratic 
in the bosonic operators b. In contrast to the expres- 
sions of Ref. the Gaussian factor e^l"''^'^/'* is absent 
in Eqs. (|35|l . and (|39|l as it is included in the interac- 
tion potential voiq) [Eq. iflfijl ]. Therefore, from Eq. 
the CF density operator p(q) has the following bosonic 
representation 

p(q) - PT('l) + ^i('l) (40) 
= <5q,oiV; + 2z^sin(qAkrV2)fe]j+k^k. 

k 

We are now able to study the complete Hamiltonian 
(|30|l . i.e., the interacting system at z^* = 1. The bosonic 
representation of the Zeeman term Hz is given l)y ' ' 

Hz=7Y.^l,b^~hK- (41) 
q 

Because the CF density operator is quadratic in the 
bosonic operators, Hj in Eq. H3Q|I is quartic in the oper- 
ators 6, and one obtains with the help of Eq. (|40|l 

Hi = T. /fJ-o(k)[Fo^^(fc)]%in(kAprV2) 

xsin(kAqrV2)&i+qMp-k^ 
= E /0«o(k)[F„^^(fc)]%in^(kAqrV2)&t6, 

+ Y.[P^-o{k) [^^o^^(fc)]%in(kAprV2) 
X sin(k A qrV2)&UqfoJ,_k^^- (42) 

Adding Eq. (|41|l to (|42|l yields the expression for the total 
Hamiltonian of the FQHE system at i^* = 1 in terms of 
bosonic operators, 

H^=H^+HP^„ (43) 
with the free boson Hamiltonian 

Hi = -\lN; + J2w<ib% (44) 
q 

and the boson-boson interaction part 

k,p,q 

xsin(kAq/*V2)&Uq&|,_kMp- (45) 
The dispersion relation of the bosons is given by 

^q = 7 + / 0z.o(k)[Fo^^(fc)]%in2(kAqrV2). (46) 




FIG. 2: Schematic representation of the neutral excitations 
of the quantum Hall ferromagnet of CFs a.t v* — 1. The 
bosonization method takes into account only spin- wave (SW) 
excitations which involve a spin reversal inside the lowest CE- 
LL. The excitations between different CF-LL, without a spin 
reversal (charge mode, CM) and with a spin-flip (SF) are 
neglected. 

Therefore, the bosonization method for the 2DEG at 
v = 1 allows us to map the fermionic model Hamiltonian 
(|30|l for the FQHE system at i^* = 1, into an interacting 
bosonic one. The ground state of the bosonic model l|43() 
is the boson vacuum, which corresponds to the quantum 
Hall ferromagnet of CFs, with energy EpM = 
In fact, this is also the ground state of the noninteracting 
bosonic model (|41|l . and thus the ground state of the 
system does not change when the interaction between 
the electrons is introduced. 

The bosonization approacbi^ restricts the total Hamil- 
tonian (|30|l to one particular sector of all excitations of 
the quantum Hall ferromagnet of CFs, namely the one 
which corresponds to a spin reversal inside the lowest 
CF-LL: a spin- wave (SW) excitation (Fig.|21). It neglects 
excitations between different CF-LLs, both those with- 
out spin reversal (charge mode, CM) and those which 
also involve a spin- flip (SF). 

A further issue in the present bosonization scheme 
would be to provide the expression of the fermion field 
operator in terms of the bosonic ones, 6^ and &q, 
i.e., the analog of the Mattis-Mandelstam formula for 
one-dimensional fermionic systems^^ii^ Such expression 
would enable one to calculate, for instance, correlation 
functions related to the model (41). Here, however, we 
are interested only in the energy dispersion of the collec- 
tive spin excitations, which find their natural formulation 
in terms of the bosonic operators and which may be cal- 
culated without the precise expression for the fermionic 
fields. 



A. Dispersion relation of the bosons 

In this section, we evaluate the dispersion relation of 
the bosons [Eq. (|46|l ] and investigate whether the inter- 
action potential between them accounts for the forma- 
tion of a bound state of two bosons, as it was found for 
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the 2DEG at u = 1. In addition to the bare Coulomb 
potential, we also consider the modification of the in- 
teraction potential between the original fermions due to 



finite-width effects. 

For the Coulomb potential [Eq. , the integral over 
momentum in l|4(i|l may be performed analytically, 



J 



7+ / g «o(k)e-l'='-^l^/2 



7+ - c2 

elB V ^ 



a 



sin2(kAqrV2) 



r 



C" |,/*cp/4(2-c2)T ( k'*c|- 



4(2 -c2) 



(47) 



where Io(x) is the modified Bessel function of the first 
kind.-2i The function 

„5 

a(c) = 1 - 2c3 + 



%/2"=~? 

depends on the charge of the pseudo- vortex, — c^, which 
is related to the electronic filling factor u. Remember 
that, for u = l/(2s + 1) we have = 2s/(2s + 1). For 
the case of vanishing Zeeman energy (7 = 0), the results 
of Eq. l|T7|) are shown in Figs.|21(a) and (b) (solid lines) 
for J/ = 1/3 and 1/5, respectively. The energies are given 
in units of c^/cIb- 

Notice that at g = the energy of the bosons is equal 
to the Zeeman energy 7 (in agreement with Goldstone's 
theorem) and that Wq ~ |(7?bP in the small- wavevector 
limit, Iqlnl *C 1 [see Eq. This is exactly the behav- 

ior of spin-wave excitations in a ferromagnetic system. 
Therefore, the long-wavelength limit of corresponds 
to spin-wave excitations of the quantum Hall ferromag- 
net of CFs. To the best of our knowledge, this is the first 
analytical calculation of Wq. 

In analogy with the quantum Hall ferromagnet a,t v = 
1 , we may associate the limit g — > 00 of Wq to the energy 
Eph to create a well-separated CF-quasiparticle-quasihole 
pair, i.e. 



Eph = 7- 



TT / 

2^ 



c^a{c)—— 



1 + 



where 



0.163 v 
0.051 1^ 



1/3, 
1/5, 



(48) 



(49) 



in units of e^/els- 

Our result is in agreement with the one found by 
Murthy in the "naive" approximation, where the dis- 
persion relation of spin-wave excitations of the FQHE 
system (neglecting disorder, finite-width and Landau- 
level- mixing effects) at J' = 1/3 has been evaluated using 
the time-dependent Hartree-Fock approximation. It is 



not a surprise that both calculations are in agreement 
as the main approximation of our bosonization method 
[Ea. H34|) ] resembles the one considered in the random- 
phase approximation (RPA)iS 

On the other hand, our results appear to be overes- 
timated when compared to the ones of Nakajima and 
Aoki, who have also calculated the spin wave excitation 
of the FQHE system at = 1/3 and 1/5.^^ They started 
from the results of Kallin and Halperin for the 2DEG 
a.t v = 1,^5 turned to the spherical-geometry approach, 
introduced the CF picture and then performed mean- 
field calculations. Their results agree with the ones of 
Ref. (33 for the 2DEG at z/ = 1 with the replacements 
e ^ e* and Ib ^ I*- Finally, Mandal has studied spin- 
wave excitations in the framework of a fermionic Chern- 
Simons gauge theory^S For 1/ = l/(2s-|-l), his results co- 
incide with those obtained by Kallin and Halperin for the 
2DEG at ly = 1, with the replacement Ib — > I*, wheras 
the charge remains unchanged. Mandal's estimates for 
the energy of a well separated CF quasiparticle-quasihole 
pair are therefore larger in comparison with the above 
mentionned works. A possible reason for this mismatch 
may be the presence of higher CF-LLs, which are im- 
plicitely taken into account in the numerical approach by 
Nakajima and AokiSi and explicitely in Murthy's time- 
dependent Hartree-Fock calculations, beyond the naive 
approximation,^- Remember that the model (|3Uf) used 
in our analyses, however, is restricted to a single level 
and thus ignores the presence of higher CF-LLs. The 
latter may in principle be included in a screened inter- 
action potential, the g-dependent dielectric function of 
which may be calculated within the RPA.^^ However, 
this would only take into account virtual density-density 
fiuctuations, but no single-particle excitations. 

The energy E^'f^ has been calculated numerically by 
Mandal and Jain for the FQHE system (neglecting dis- 
order, finite width and Landau level mixing) dA, v — 
1/3, using trial wavefunctions for the ground state in 
the spherical-geometry approach. "^^ The thermodynamic 
limit obtained from the calculations on finite systems in- 
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limit, and the inclusion of finite- width effects cuts off the 
short-range contributions of the Coulomb potential. 

From the experimental point of view, spin reversal ex- 
citations of the FQHE system have been studied with the 
help of inelastic light scattering. In this case, it is 
possible to probe some critical points of the dispersion re- 
lation of the neutral excitation, namely the ones at 5 — > 
and 5 — > cxD. In particular, the latter can be probed be- 
cause the wavevector is no longer conserved in the pres- 
ence of residual disorder. The experimental spectrum 
at i/ = 1/3 consists of five different peaks. ''^■^^ Three of 
them have been associated to the charge modes at q ^ 
(Ao), q ^ CO (Aoo), and to the magnetoroton minimum. 
The one at the Zeeman energy has been associated to the 
long-wavelength limit of the spin- wave excitations, while 
the one between Aq and Aqq was related to Eph- In this 
case, it was estimated that i?^ = 0.054 e'^/els- We ob- 
tain precisely this value for a width parameter A = 3.0, 
which corresponds to 290A at B = 7T {Is = 97 A). In 
view of the approximations made (no Zeeman splitting, 
restriction to a single CF-LL level, and neglecting impu- 
rity effects) , this is surprisingly close to the experimental 
situation, where a quantum well of width 330 A was in- 
vestigated at 7T4i*^ 



FIG. 3: Dispersion relations of the bosons as a function of 
the momentum q with 7 = 0, in units of e^/els, when (a) 
v = 1/3 and (b) v = 1/5. The solid line corresponds to Eq. 
H47|l . whereas the dashed and solid gray lines represent the 
results for the width parameters A = 0.5 and 1.0, respectively. 

dicated that — 0.074 e^/elB, in agreement with the 
previous results by RezayijS who performed numerical- 
diagonalization calculations on finite systems, without 
the CF picture. Notice that this value is almost half 
as large as the one we obtain in the framework of the 
Hamiltonian theory [Eq. (|^ ]. 

In order to include finite-width effects, we replace the 
bare Coulomb potential (|16|l by 

= ^e-\'^iB\V2^\gM' [1 _ Erf(Ag)] , 

where Erf(a;) is the error function. '^^ This expression has 
been obtained under the assumption that the confining 
potential in the z-direction (with a characteristic width A, 
in units oUb) is quadraticiSS*^ After replacing Vo{q) 
VQ{q), the momentum integral in Eq. (|46ll is evaluated 
numerically. For v = 1/3 and 1/5, the results for two 
different values of the width parameter (A = 0.5 and 1.0, 
dashed and gray lines, respectively) are shown in Figs. |3| 
(a) and (b). 

Notice that, in agreement with results by Murthy, the 
energy of the bosons decreases as A increases*^ In fact, 
the finite-width results are likely to be more reliable than 
the ones at A = as the preferred combination for the 
charge density [Eq. (|26|l ] is derived in the long-wavelength 



B. Bound States of Two Bosons 

Apart from the fact that our bosonization scheme re- 
produces the time-dependent Hartree-Fock (or RPA) re- 
sults, it goes beyond as it yields an interacting bosonic 
model (|43f) . Such kind of result was first derived 
for the 2DEG at 1/ = 1 in Ref. ITsl In this sec- 
tion, we follow the lines of this work and study two- 
boson states in the interacting bosonic model in order 
to find out a relation between the two-boson bound 
states and the quasiparticle-quasihole (small skyrmion- 
antiskyrmion) pair excitations. 

Finite-width effects are not considered here, and we 
therefore study the bosonic model H43() with the disper- 
sion relation of the bosons given by Eq. H47|) . In or- 
der to study two-boson states, one needs to solve the 
Schrodinger equation 

H|$p) = ^p|$p), (50) 

where |<i>p) is the more general representation of a two- 
boson state with total momentum P, namely 

I'&p) = E*P(q)^f-q4+J^Af). (51) 

q 

Here, P is a good quantum number because the total 
momentum of a 2D system of charged particles in an 
external magnetic field is conserved when the total charge 
of the two-particle state is zero.^^ 

Substituting Eq. H51|l into H50|) we obtain the following 
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eigenvalue equation 

[e-£;p(q)]$p(q) = yd2fcAp(k-q,q)$p(k), (52) 

where £'p(q) = wp/2-q + iL'p/2+q is the energy of two 
free bosons, e = Ep—EpM, and the kernel of the integral 
equation reads 

In the above expressions, all momenta are measured in 
units of the inverse of the magnetic length l'^^ . 

The eigenvalues of H52|l are evaluated numerically, us- 
ing the quadrature technique. ■^^ In this case, the integral 
over momenta in H52|) is replaced by a set of algebraic 
equations, 

[s - £;p(qi)]$p(qi) w J2 Cj-^p(qj - qi,qi)*p(qj), 

(54) 

and the eigenvalues for a fixed momentum P may be 
calculated by means of usual matrix techniques. The co- 
efficients of the quadrature Cj and the points q^ depend 
on the parametrization. We chose the one with 61 points, 
which is used to calculate two-dimensional integrals over 
a circular region.^ A cut-off for large momenta has to 
be introduced in order to define the region of integra- 
tion. Therefore, only bosons with {qlsl < 1 are included 
in our calculations, and we concentrate on the analysis 
of the lowest-energy two-boson state. Both restrictions 
are related to the limitations of the numerical method. 

In Figs. 21 (a) and (b), we show the lowest-energy state 
of two bosons (solid line) as a function of the total mo- 
mentum P for v = 1/3 and 1/5, respectively, in units 
of e^/els, in the limit of 7 = 0. The shaded area corre- 
sponds to the continuum of scattering states. Because the 
solid line is below the band of scattering states, we may 
interpret the lowest energy eigenvalue (Ep) as a bound 
state of two bosons. We have also included (dashed 
line) the value of the energy E^^ of a CF-quasiparticle- 
quasihole pair [Eq. H48|) ]. with the particle and the hole 
far apart from each other. Notice that in both cases, the 
asymptotic limit (|P| — > 00) of Ep is close to i?^. In 
analogy with the 2DEG at i' = 1,^^ the bound state of 
two bosons with large P may be understood, for instance, 
as a CF quasiparticle bound to one spin wave plus a dis- 
tant CF quasihole.^'^ In this scenario, our results indicate 
that the creation of a CF quasiparticle dressed by a spin 
wave (or alternatively, a small skyrmion of CFs) is more 
favorable for the FQHE system at = 1/3 than for the 
one at z/ = 1/5 as |£'^j — E'p^ool is smaller in the lat- 
ter case. At present, no Knight-shift data are available 




FIG. 4: Dispersion relations of the two-boson state as a func- 
tion of the total momentum P, with 7 = 0, in units of e^/eZs, 
when (a) v = 1/3 and (b) 1/5. The solid line corresponds 
to the lowest-energy state of two bosons, the dashed line is 
the energy of a well-separated CF quasiparticle-quasihole [Eq. 
II48II 1. and the shaded area represents the continuum of scat- 
tering states. 

for FQHE systems at z^ = 1/5 to check our theoretical 
results. 

The energy of small CF skyrmions for FQHE systems 
at V = 1/3 and 1/5 has been calculated numerically by 
Kamilla et al.^ using a hard-core wavefunction4& How- 
ever, it is not possible to compare our results, namely, the 
asymptotic limit of Ep , to the ones obtained by Kamilla 
and co-workers as they have evaluated the energy only 
of the positively charged excitation. In contrast to the 
2DEG aX V = 1, the particle-hole symmetry is broken 
ai V — 1/3 and 1/5. Therefore, the energy of the posi- 
tive excitation cannot be used to estimate the energy of 
the negative one, and consequently the energy of a well- 
separated skyrmion- antiskyrmion pair cannot be deter- 
mined. The breaking of the particle-hole symmetry at 
V = l/i may be observed in the Knight-shift data, which 
are asymmetric around this filling factor^ 



IV. SKYRMIONS OF CFS 

It has been shown in Rcf. 13 that the semi-classical 
limit of the interacting bosonic model provides a micro- 
scopic basis for a phenomenological model used by Sondhi 
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et ali^ for the description of quantum Hall skyrmions. 
Sondhi's model is given by the Lagrangian density 



^eff = ^PoA{n) ■ dtn - ^ps{Vnf 



1 



:7Po/^sn • B 



(55) 



where ps is the spin stiffness?*? A{n) is the vector po- 
tential of a unit monopole (e°'"^9a-4h = n'^) and pa — 
i'/{2TTlg) is the average electronic density. The topolog- 
ical charge or skyrmion density is given by 



q{r) ^ ^e°''^n ■ {daii x d^n), 
an 



(56) 



with a,b,c = x, y, z and a, /3 = cc, y. Within the 
bosonization approach, the gradient, Zeeman and in- 
teraction terms of the above Lagrangian density have 
been derived, and one obtains, from a microscopic de- 
scription, exactly the same coefficients. This result cor- 
roborates the relation between the bound state of two 
bosons and the small skyrmion-antiskyrmion pair excita- 
tion discussed in Ref. il^ In other words, the interacting 
bosonic model may be used to describe the quantum Hall 
skyrmion at i/ = 1. 

Here, we proceed in the same manner in order to derive 
the analogue of Sondhi's model for the CF skyrmion at 
V = l/(2s-|- 1). We analyze the bosonic model H43|l . with 
the dispersion relation of the bosons given by Eq. (|47|l . 

Following the lines of Ref. |Oj we describe the CF 
skyrmion by the coherent state 



sk'-'^ (n)) = exp(-iO)lFM), 



(57) 



where the operator O is a non-uniform spin rotation 
which reorients the local spin from the direction to 
n(r), 



O 



■ S(r) 



d^q [l]+(q)5rq + f7-(q)5j:^ 



(58) 



Here, S(r) is the spin operator, n(r) is a unit vector, 
and r2(r) = x n(r) defines the rotation angle. We 
assume that r2(r) describes small tilts away from the 
direction, i.e. ri*(q) vanishes when \qlB \ ^ 1- 

Using the expression for the spin density operator in 
terms of the boson operator the operator O becomes 



q 



(59) 



and the constant /3 will be defined later. 

We start with the analysis of the CF density operator 
(|40|l . The average value of this operator in the state 1)5 7|l 
is 

X ^ sin (k A qrV2) f^k+qf^^q- (60) 



Because only the long-wavelength limit of the theory is 
considered, one may approximate 

sin (k A qrV2) « k A qrV2 = (k x q)^/*V2, 



and the CF form factor F^p{k) w 1 — c^. 

This implies that we only retain terms 0{k^) in Eq. 
(|60|l . Replacing the sum over momenta by an integral, 
and using the relation between the vector r2(r) and the 
unit vector field n(r), the average value of the CF density 
may be written as 



(k + q)r!k+q X qr!+ 



= -(1 - c2)27riV;/32(Z*)4 j (fr e-'^ ^'e, ■ (Vn^ x Vn^' - x Vn^) 



and therefore, its Fourier transform becomes 

(^(r)) = -(1 - c^)2TTNl(3^{r fe, ■ (Vn=^ x Vn^ - Vn" x Vn^) . 

I 



(61) 



(62) 



Notice that, if we choose the constant [3 — is equal to the topological charge times the filling factor 

1/ Uir.flfArf) the average value of the CF density = (1 ^) = V(2s + 1), as suggested by Sondhi et al. 
^ V y in Ref. llOl With this choice, we may calculate the aver- 
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age value of the energy in the state H57|) , considering the 
interacting bosonic model (|43|) . with Wq as in Eq. 147|) . 

The average value of the quadratic term of the Hamil- 
tonian (|43() is given by 



Because the long-wavelength limit of the dispersion rela- 
tion of the bosons [Eq. H47|) ] is 



~ 7 + 1 Vl 



e 



TT 1 



2 4 



a{c)\ql 



* |2 



with the constant 



a(c) = 1 - 2c-'^ + 



(2-c2)3/2^ 
the average value of Hq is 

(Wo) « ips / [Vn(r)]^ 



(63) 



(64) 



Here, the spin stiffness is defined as 



1 /tt 



^^=16^d^V2"(^^^^^''^'^' 



with 



Ps 



1 e*^ hn 



IQtt el* V 2 



(65) 



(66) 



(67) 



and /i(c) = a(c)/(l — c^)^. Notice that one would have 
obtained the value p*g for the CF spin stiffness simply by 
replacing e ^ e* and Ib — > I* in the expression for the 
electronic spin stiffness ai v = This corresponds to 
a CF mean-field approach, ■^^ and /i(c) may thus be inter- 
preted as a mismatch factor with respect to the mean- 
field result. 

Performing an analogous calculation for the quartic 
term of the Hamiltonian H43|l , one finds 

{Hi) ^ MJ2 Mk) [Fo^^(fc)]%in(kAprV2) 

k,p,q 

X sin(k A qrV2)r!^+qr!p_^r!tpOtq 



d^rd^r' V{r-r')q{r)q{r'), 



(68) 



where TV = 2/ (A{'i7r)'^Nf{l*)^y Therefore, using 

equations H65|l and (|68|l . one obtains the average value 
of the total Hamiltonian with respect to the state Isfc*-^^) 

{H) = I d\ [Vn{v)f + ^g*pBB^N; 



+ 



e2(l-c2)2 



(69) 



Eq. (|69|l corresponds to an energy functional of a static 
configuration of the vector field n(r), derived from a La- 
grangian density similar to (|55|l . Notice that the coeffi- 
cient of the Zeeman term. 



9 fJ-B 

2 27r/*2 



9 IJ-B V 
2 27r?| ' 



is identical to that of Sondhi et al.,^'^ with the differ- 
ence that here, this term is derived from first principles 
and not assumed phenomenologically. Furthermore, the 
last term in Eq. (|59|l accounts for the Coulomb inter- 
action between CF-quasiparticles with effective charge 
e* =e(l -c2). 

The analysis of the semi-classical limit of the interact- 
ing bosonic model H43I) allows us to calculate analytically 
the spin stiffness. From Eq. H66I) . one obtains 



PS 



6.207 x 10-3 iy=l/3, 
2.269 x 10-3 iy=l/5. 



(70) 



in units of e^/elB- 

There is a set of different values for the spin stiffness at 
1/ = 1/3 and 1 /5 available in the literature. Based on Ref. 

Mandal and Ravishankai4^ estimated that the spin 
stiffnesses at = 1/3 and 1/5 are equal to 1.591 x 10"^ 
and 4.774 x 10—* e'^/els, respectively. These values are 
derived from the spin stiffness ad, v — 1 with the replace- 
ments e e* and Ib — > I*, and may be interpreted 
as the mean-field values of the CF spin stiffness, as in- 
dicated above. In fact, these results are in agreement 
with estimates by Sondhi et ali^ The spin stiffness has 
also been calculated numerically by Moon et ai.jii us- 
ing the hypernetted-chain-approximationjS& In this case, 
it was found that ps = 9.23 x 10"'* [v = 1/3] and 
2.34 X lO-'' eV(e^B) W = 1/5]- Those values are also 
considered by Lejnell et al.^^ in the study of the quan- 
tum Hall skyrmion. Furthermore, the spin stiffness at 
ly — 1/3 as a function of temperature and finite width has 
been estimated by Murthy^^ Hartree-Fock approx- 

imation of the Hamiltonian theory. At low temperatures 
{0<T <2K), one finds psiT) « 2.125 x lO^^. 

The values we obtain with the help of the bosoniza- 
tion approach in the Hamiltonian theory are larger than 
the theoretical estimates found in the literature. In the 
case of the mean-field values obtained by Sondhi et ali^ 
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and Mandal and Ravishankar,— the mismatch is given 
by the factor ^(c), which is 3.9 a.i v = 1/3 and 5.1 at 
V — 1/5. The origin of this mismatch is the Gaussian 
exp(— |(7/bP/2c^) in the CF form factor l(77|) . which takes 
into account the difference between the magnetic length 
of the electron and that of the pseudo- vortex, which con- 
stitute the CFi The mismatch is thus due to the internal 
structure of the CF, which is not taken into account in 
a mean-field approximation, where the CF is considered 
as a pointlike particle with charge e*. 

A comparison with experimental data has only been 
done in Ref. ls^ where the temperature dependence of the 
spin polarization at u = 1/3 has been evaluated, using 
the continuum quantum ferromagnet model in the large- 
N approach. In this work, the two free parameters of 
the model (spin stiffness and magnetization density) have 
been taken temperature-dependent. This formalism had 
previously been used to describe the spin polarization 
of the 2DEG sX v = 1^ but with constant parameters. 
The agreement between the theoretical result and the ex- 
perimental data^ even when only the zero-temperature 
values of two free parameters were taken into account, in- 
dicated that, within this formalism, the estimated value 
of the spin stiffness is quite reasonable. Although we 
find a larger value for the spin stiffness than Murthy's, it 
can be reduced if finite-width effects are included in the 
model. Remember that the preferred combination of the 
CF density is less accurate at larger momenta and the 
inclusion of finite-width effects softens the short-range 
terms of the interaction potential. 

Finally, the energy of the CF skyrmion {Esk) may also 
be estimated in the long- wavelength limit. The leading 
term of the skyrmion energy is Ai:ps which is basically 
the energy of the skyrmion in the limit of zero Zeeman 
energy. One obtains from the values (|70|l 



Esk — 



0.078 
0.028 



1/3, 
1/5, 



(71) 



in imits of /els- These results are in agreement with 
the energy of a large (number of spin reversals) skyrmion 
calculated by Kamilla et alJ^^ This indicates that the 
hard-core wavefunction description for the CF skyrmion 
and the Hamiltonian theory in the bosonization approach 
converge in the long- wavelength limit. 



V. SUMMARY 

In conclusion, we have calculated neutral spin exci- 
tations of FQHE systems at i/ = l/(2s + 1), using a 
bosonization approach for the Hamiltonian theory of the 
FQHE proposed by Murthy and Shankar.-''. The general- 
ization of this formalism with the inclusion of a discrete 
degree of freedom has been presented, under the assump- 
tion that one may neglect higher CF levels, which are 
separated from the ground state by a larger gap than the 
spin-level splitting (7 — > limit). 



Because in the CF picture the electronic system at 1/ = 
1 / (2s -I- 1) is described by a CF model with effective filling 
factor ly* — 1, we can apply the bosonization method for 
the 2DEG at = 1, recently developed by two of us.~ In 
this case, the neutral spin excitations of the system are 
treated as bosons, and the original interacting CF model 
is mapped onto an interacting bosonic one. 

The dispersion relation of the bosons (neutral spin re- 
versal excitations) has been evaluated analytically for 
an ideal system, neglecting finite-width. Landau-level- 
mixing and disorder effects. We have illustrated the re- 
sults for FQHE systems at ly — 1/3 and ly = 1/5. The 
bosonization approach allows us to calculate analytically 
the energy of the spin excitations over the whole momen- 
tum range. 

Our results agree with numerical investigations in the 
time-dependent Hartree-Fock approximation.^- We have 
shown that the long-wavelength limit of the bosonic dis- 
persion relation corresponds to spin-wave excitations of 
the quantum Hall ferromagnet of CFs, which is the 
ground state of the system. The energies of the spin 
waves are larger than those obtained in previous the- 
oretical studies by Nakajima and Aokii^ The large- 
wavevector limit of the dispersion relation is related to 
the energy to create a CF-quasiparticle-quasihole pair 
{Eph), the constituents of which are far apart from each 
other. Our results are also larger than numerical ones 
obtained with the help of trial wavefunctions for the ex- 
cited state in the CF picture^ and exact-diagonalization 
studies on systems with a few number of particlesJS 
Even if these values occur to be closer than ours to ex- 
perimental estimateg^^'^^ of the Eph, an excellent agree- 
ment is obtained when finite-width effects are taken into 
account. Indeed, the width of the quantum well used 
in the experiments corresponds to a width parameter 
X 31b (~ 300 A at 7T), for which we obtain a reduc- 
tion of Epfi roughly by a factor of three. Note that finite- 
width results are more reliable in the Hamiltonian the- 
ory, which was originally derived in the long-wavelength 
limit, than that for the ideal 2D case because, in this case, 
short-range contributions of the interaction potential are 
suppressed 1^ This cuts off short- wavelength fluctuations 
and thus brings the model back to its regime of validity. 

In the framework of the interacting bosonic model, we 
have shown that the interaction potential between the 
bosons accounts for the formation of two-boson bound 
states. As in the case of the 2DEG at v = 1, the two- 
boson bound state is interpreted in terms of a small CF- 
skyrmion-antiskyrmion pair excitation. FQHE systems 
at 1/ = 1/3 and 1/5 have been investigated. Based on the 
relation between the asymptotic limit of the energy of 
the lowest two-boson bound state and Eph, we conclude 
that a CF quasiparticle dressed by spin waves is more 
stable than a bare CF quasiparticle only at v = 1/3. 
Experimental results at v = 1/5, such as Knight shift 
data, are lacking at the moment and may give further 
insight into the spin-excitation spectrum of this state. 

Finally, we have shown that the semi-classical limit of 



14 



the interacting bosonic model yields an energy functional 
derived from a Langrangian density, similar to the phe- 
nomenological one considered by Sondhi et alm^ for the 
study of the quantum Hall skyrmion. We thus corrobo- 
rate the relation between the two-boson bound state and 
small skyrmion-antiskyrmion pair excitations. 
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